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ABSTRACT 


The  long-range  objective  of  the  research  is  to  develop  the  capability 
of  assessing  the  def lagration-to-detonation  transition  (DDT)  hazard  in 
HMX-based  propellants.  The  approach  is  based  on  the  concept  that  a basic 
understanding  of  the  physical  and  chemical  processes  involved  in  DDT  is 
necessary  to  achieve  this  objective.  Accordingly,  theoretical  and  exper- 
imental studies  were  undertaken  to  elucidate  mechanisms  of  DDT,  to  estab- 
lish conditions  for  its  occurrence,  and  to  formulate  a satisfactory  model 
for  its  quantitative  description. 

The  theoretical  studies  were  concerned  with  the  thermo-hydrodynamic 
treatment  of  DDT.  The  central  problem  is  to  model  the  propagation  of 
the  unsteay  flame  and  the  flow  it  produces  ahead  of  the  flame  front.  The 
flame  was  treated  as  a reactive  discontinuity  and  flow  conditions  for 
propagation  of  steady  and  unsteady  flames  were  established.  Equations 
for  the  initial  flow  produced  by  the  onset  of  deflagration  in  a closed 
tube  were  derived  to  demonstrate  that  an  accelerating  flame  produces  a 
compression  wave  ahead  of  the  flame  front.  The  initial  attempt  to  con- 
struct an  explicit  solution  for  an  accelerating  flame  in  a closed  tube, 
based  on  the  group  of  transformations  admitted  by  the  flow  equations,  was 
unsuccessful . 

The  experimental  studies  were  delayed  by  problems  of  finding  a source 
of  propellant  samples  suitable  for  Lagrange  gage  experiments.  A source 
of  HMX-based  propellant  was  found  at  Edwards  Air  Force  Base,  and  experi- 
mental procedures  were  changed  to  eliminate  the  machining  and  grooving 
operations  usually  employed  in  constructing  assemblies  for  Lagrange  gage 
experiments.  A target  assembly  was  designed  to  allow  incorporation  of 
the  gages  and  ionization  pins  into  the  propellant  charges  during  the 
casting  process. 
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I INTRODUCTION 
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Research  on  the  def lagration-to-detonation  transition  (DDT)  is 
important  to  the  Air  Force  because  many  of  the  energetic  propellants 
required  for  present  and  future  long-range  delivery  systems  are  explo- 
sively filled  compositions  that  are  capable  of  undergoing  DDT.  HMX- 
based  propellant,  for  example,  usually  bums  reliably  in  rocket  motors 
but  has  detonated  and  destroyed  the  motor  on  several  occasions.  Possible 
steps  leading  to  detonation  in  the  rocket  motor  are  fracture  of  propel- 
lant ahead  of  the  flame  and  the  subsequent  formation  of  shock  waves 
produced  by  the  increased  burning  rate  of  fractured  propellant.  But 
even  in  this  case,  the  mecahnism  of  DDT  and  the  conditions  for  its 
initiation  are  not  adequately  understood. 


The  long-range  objective  of  the  present  research  program  is  to 
develop  a computational  capability  for  assessing  the  DDT  hazard  in 
explosively  filled  propellants.  The  program  is  based  on  the  concept 
that  a basic  understanding  of  the  physical  and  chemical  processes 
involved  in  DDT  is  necessary  to  achieve  this  objective.  Combined  theo- 
retical and  experimental  studies  to  determine  the  pressure  fields  behind 
and  ahead  of  the  flame  and  to  establish  mechanisms  of  DDT  are  required 
to  develop  such  understanding. 
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II  THEORETICAL  STUDIES 

The  first  task  in  the  theoretical  study  is  to  develop  a basic 
physical  understanding  of  the  thermo-hydrodynamic  processes  governing 
the  flame  and  its  propagation.  Specific  steps  undertaken  to  attain  such 
an  understanding  are  as  follows: 

(1)  Choose  a model  for  the  flame 

(2)  Derive  and  physically  interpret  the  equations  governing 
flame  propagation 

(3)  Construct  solutions  for  the  unsteady  flow  produced  by  a 
flame  propagating  in  a closed  tube. 

It  is  convenient  first  to  choose  the  simplest  model  and  treat  the  flame 
as  a reactive  discontinuity.  Implicit  in  this  choice  are  the  assumptions 
that  the  burning  reaction  is  fast  and  that  the  hydrodynamic  flow  produced 
by  an  unsteady  flame  is  not  appreciably  influenced  by  the  detailed 
structure  of  the  flame.  In  this  case,  states  connected  by  the  burning 
process  are  governed  by  the  Rankine-Hugoniot  jump  conditions  expressing 
the  conservation  of  mass,  momentum,  and  energy  across  the  reactive 
discontinuity. 

Let  v denote  specific  volume;  F,  flame  velocity;  u,  particle 
velocity;  p,  pressure;  c,  sound  speed;  and  h,  specific  enthalpy;  and 
let  the  subscripts  i and  f denote  the  initial  and  final  states,  respec- 
tively connected  by  the  flame.  The  states  across  the  discontinuity  are 
related  by  the  Rankine-Hugoniot  jump  conditions 

vf (F-u^)  = v1(F-uf) 

(uf  " u1)2  = (pf  - p4)  (v±  - vf) 

2 (hf  - hjL)  = (pf  - Pi)  (vf  + vt) 


2 

^ . . . Ai 


(1) 

(2) 

(3) 


Equation  (1)  expresses  the  conservation  of  mass,  Eq.  (2)  expresses  the 
conservation  of  mass  and  momentum,  and  Eq.  (3)  expresses  the  conservation 
of  mass,  momentum,  and  energy.  The  locus  of  states  attainable  from  a 
given  initial  state  (p{,  v^,  h^)  lie  on  a curve  in  the  (p,v)  plane 
obtained  by  eliminating  h between  Eq.  (3)  and  the  h *=  h(p,v)  equation  of 
state  of  the  burnt  products.  This  curve,  called  the  Hugoniot  curve 
centered  on  (p^>  v^),  has  two  branches  when  the  discontinuity  is  reactive. 

When  the  reaction  is  exothermic,  the  branch  where  p.  > p.  and 

’ f i 

vr  < v.  is  called  the  detonation  branch,  and  the  branch  where  p,.  < p. 

f = i > rf  = ri 

and  Vj.  > v^  is  called  the  deflagration  branch.  With  the  present  model 
of  the  flame,  burnt  states  are  represented  by  points  on  the  deflagration 
branch  of  the  Hugoniot  curve.  Moreover,  the  maximum  flame  speed  is 
defined  by  the  Chapman-Jouguet  (CJ)  point  where  F = u^  + c^,  the  flow 
is  sonic,  and  the  Rayleigh  line  through  the  initial  state  (p  ,v^). 


2 


is  tangent  to  the  Hugoniot  curve. 

Equations  for  deflagration  CJ  parameters  will  be  derived  and 
then  the  differential  equation  governing  the  propagation  of  a reactive 
discontinuity  will  be  formulated.  Equation  of  state  information  must 
be  specified  to  compute  CJ  parameters  and  states  on  the  deflagration 
Hugoniot  curve.  Let  the  subscript  o denote  the  standard  state  and  let 
the  superscripts  x and  p denote  propellant  and  propellant  products, 
respectively.  The  initial  enthalpy  of  the  propellant  and  the  h = hP(p,v) 
equation  of  state  of  the  products  are  written  as 

h*  = AH*  + h*  (pi,v1)  (5) 

and  hP  = HP  + (6) 

o k-1 

where  HP  = AHP  - k p v /(k-1),  AH  denotes  the  heat  of  formation,  and 
o o o o o 

k denotes  the  polytropic  index  of  propellant  products.  The  combination 
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of  Kqs.  (3)  and  (6)  gives  Jgie  equation  for  the  Hugoniot  curve  centered 
on  (Pjv1)  as  * 

pf(Uvf  - v^  + Pi(vi  + vf)  = 2 q (7) 

where  y = (k  + 1)/  (k  - 1)  and  q = [ (AH*  - AHP)  + hX(p^,v^)  + k p^v^/k-]]. 

At  this  stage  it  is  convenient  to  neglect  the  kp  v /(k-1)  term  in 

o o 

q and  introduce  the  nondimensional  variables  pf/p.t.=  P»  v /v  = v» 

2 2 i i i j- 

q/p^v^  = ^ and  U = ^Uf  ~ Ui^  ^PiVi‘  Ecluations  (2)  and  (7)  can  then  be 
written  as 


(P-1)  (1-V)  = U2  -(8) 

and 

P(yV-l)  + (V+l)  = 2Q  (9) 

2 

Using  Eq.  (4)  and  the  definition  of  sound  speed  c = kpv, 
tion  F-u^.  = c^  can  be  written  as 

£=i  = k I 

l-V  V 

'Hie  CJ  parameters  are  obtained  by  solving  F.qs.  (8)  - (10) 

The  combination  of  Eqs.  (8)  - (10)  leads  to  the  following 
relating  U and  V to  Q. 

pU2  + (p  + 1)  - 2Q  - 0 (13) 

(1-V)2  + 2B(1-V)  - (H  -)  B = 0 (12) 


the  CJ  condi- 


(10) 

and  Eq.  (1). 
equations 


where  B = (2Q/(y+l)  - 1).  The  roots  of  Eqs.  (11)  and  (12)  give  the 
particle  velocities  and  volumes  at  the  CJ  points  on  the  detonation  and 
deflagration  branches  of  the  Hugoniot  curve  centered  at  (P^j)*  Solving 
Eq.  (11)  and  converting  back  to  dimensional  variables  gives  the  equation 


uf  = ui  ± k^T  [2{(k2-l)q  - k(k+l)p1v1)}]1/2 


(13) 
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where  the  positive  sign  gives  the  particle  velocity  at  the  CJ  point  on 
the  detonation  branch  of  the  Hugoniot  and  the  negative  sign  gives  the 
particle  velocity  at  the  CJ  point  on  the  deflagration  branch. 

Whereas  liberation  of  chemical  energy  in  a detonation  produces  mass 
motion  directed  toward  the  reaction  front,  the  liberation  of  chemica] 
energy  in  a deflagration  produces  mass  motion  directed  away  from  the 
reaction  front.  Expanding  the  square  root  term  in  the  equation  for  the 
roots  of  Eq.  (12)  gives  the  following  equations  for  the  CJ  volumes  on 
the  detonation  (14)  and  deflagration  (15)  branches  as 


and 


v ~ k+1 

If.  = 2(k-l)q  _ J<__ 

v^  k k+1 


(14) 


(15) 


Rearranging  Eq.  (1)  to  give  the  equation 


F - u. 


<Vui> 

(V-l) 


(16) 


2 2 

and  eliminating  (u^-u^)  by  using  Eq.  (13)  with  A = 2((  k -l)q  - 
k(k+l)p^v^)  gives  the  equation  for  the  maximum  deflagration  velocity  as 


F " ui  = 2 

(1+A  /kp^) 


(17) 


The  equation  for  the  CJ  deflagration  pressure  then  follows  from  Eq.  (9) 


P-1 


A /Plvt 


(k+1)  (1+A/Ap1v1) 


(18) 


when  A'VkpjVj  »1,  Eq.  (17)  gives  the  maximum  deflagration  velocity  as 

(19) 


PiVi 

F - u = k 

i A 


and  the  corresponding  equation  for  the  change  in  pressure  across  the 
discontinuity  follows  from  Eq . (19)  as 

kpl 

(20) 


Pi  " Pf  “ k+T 
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Equations  (13),  (15),  and  (18)  show  that  the  reaction  in  a CJ 
deflagration  wave  produces  a decrease  in  particle  velocity  and  pressure 
but  an  increase  in  volume.  Other  states  on  the  deflagration  Hugoniot 
curve  must  be  considered,  however,  because  the  CJ  deflagration  represents 
a hypothetical  idealized  case.  These  states  lie  on  the  weak  deflagration 
branch  where  p^,j  < — P^>  an<^  where,  according  to  Jouguet's  rules,  the 

flow  of  reaction  products  is  subsonic  with  respect  to  the  reaction  front. 
Jouguet's  rules  for  the  flow  in  front  of  a weak  deflagration  give  an 
insight  into  the  problem  of  modeling  unsteady  deflagration  and  DDT.  The 
flow  ahead  of  the  wave  is  influenced  by  the  wave  itself  because  the  flow 
ahead  of  the  wave  is  subsonic  with  respect  to  the  wave  front.  Consequently, 
an  accelerating  flame  produces  an  unsteady  compression  wave  ahead  of 
itself  and  the  initial  pressure  for  the  burning  process  increases  as  the 
flow  develops.  The  burnt  states  are  represented  by  points  on  the  family 
of  deflagration  Hugoniot  curves  centered  on  an  adiabatic  compression 
curve  of  unburnt  material.  Moreover,  the  continued  acceleration  of  the 
flame  will  result  in  shock  formation  and  initiation  of  detonation  in 
material  ahead  of  the  flame  front. 


The  differential  equation  governing  the  propagation  of  a reactive 
discontinuity  will  now  be  derived  to  obtain  a more  quantitative  description 
of  unsteady  deflagration  waves.  Let  t and  x denote  time  and  Eulerian 
distance,  respectively,  and  p^  - p^.  = Ap  and  u.  - u^.  = Au.  When  the 
deflagration  propagates  as  a discontinuity,  the  time  variations  of  p^ 
and  Pj.  satisfy  the  equations 


= 3e  + F 3r 

Dt  3t^  Sx^ 


(21) 


!!!l  . M + F 9p 

Dt  3tj  9Xj. 


(22) 


where  the  subscripts  i and  f denote  that  the  partial  derivatives  are 
evaluated  at  the  top  and  bottom  of  the  discontinuity.  Subtracting 
Eqs.  (21)  and  (22)  after  the  partial  time  derivatives  have  been  eliminated 
using  the  identity 
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. . ^P  + u IP 
dt  at  ax 


(23) 


leads  cj  the  equation 


5£P  = d£ 

Dt  dt 


+ (F-u1) 


(F_Uf)^ 


(24) 


A similar  procedure  leads  to  the  corresponding  equation  for  u 


DAu  du  du  . . au 

Dt~  = dl.  " dtf  + (F-Ui)a^1 


s du 
(F-u, ) 

f dx. 


(25) 


The  equation  governing  the  discontinuity  is  obtained  by  combining  Eqs. 
(24)  and  (25)  with  the  equations  of  motion,  which  are  written  as 


n—  = _ IP 

ydt  ax 


(26) 


d£  _ 

dt 


2 

c p 


3 n 
dx 


(27) 


Multiplying  Eq . (25)  by  m = p^(F-u^)  = p^(F-Uj-)  and  adding  the  resulting 
equation  to  Eq . (24)  leads  to  the  equation 


Mp  + 

Dt 


DAu 

Dt 


~ P1(cj  -(F~ui)2) 


3u 
3x . 


£ <28> 


Equation  (28)  gives  hydrodynamic  conditions  associated  with 
deflagration  waves.  Consider  first  a deflagration  wave  propagating  at 
constant  velocity  when  D(Ap)/At  = D(Au)/Dt  = 0 and  both  terms  on  the  right- 
hand  side  of  Eq . (28)  are  zero.  Necessary  conditons  for  this  case  are 
either  du/dx^  = du/dx^  = 0,  or  du/dx^  = 0 and  cf  = F - u^,  because  unreacted 
material  satisfies  the  condition  c^  > F - u . Solutions  involving  a pre- 
compression shock  followed  by  a deflagration  wave  can  be  constructed  to 
satisfy  both  of  these  conditions.  Rejection  of  the  CJ  condition  on 
physical  grounds,  however,  gives  du/ dx^  = du/dx^.  = 0 as  the  necessary 
hydrodynamic  condition  for  steady  deflagration  waves.  Of  particular 
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interest  is  the  flow  when  the  particle  velocity  of  burnt  material  is 
zero.  In  this  case  the  strengths  of  the  shock  and  deflagration  waves 
are  such  that  the  increase  in  particle  velocity  produced  by  the  shock 
is  exactly  offset  by  the  decrease  in  particle  velocity  produced  by 
the  deflagration. 

Now  consider  the  case  of  an  accelerating  deflagration  D(Ap)/Dt  > 0, 
D(Au)/Dt  > 0 supporting  a compression  wave.  Because  c.  > F-u.  and 
du/bx^  < 0 in  the  compression  wave,  and  c^  > F-u  in  the  products,  the 
possibility  arises  that  Eq.  (28)  can  be  satisfied  when  Jju /ax^.  is  positive, 
zero,  and  negative. 

The  flow  produced  by  the  propagation  of  a deflagration  wave  from  a 
rigid  wall  is  of  particular  interest  because  our  experimental  study  of 
DDT  is  designed  to  model  this  rear-boundary  condition.  Equations  for 
the  initial  flow  will  be  derived  for  this  casa  to  show  that  a compression 
wave  must  be  formed  as  the  deflagration  propagates  from  wall.  Let  the 
superscript  o denote  the  initial  condition  at  time  t = 0.  Then  initially 
the  rear-boundary  conditon  is  expressed  by  the  equations 
Uj-  = (duj./dt)°  = 0,  We  also  assume  that  the  unburnt  material  is 
initially  at  rest,  u^°  = 0,  and  that  the  initial  velocity  of  the  flame 
is  zero,  F°  = 0.  It  follows  from  the  jump  conditions,  Eqs.  (1),  (2), 
and  (3),  the  pf°  = p^0,  and  h^°  = h^0,  but  that  pj.°  <<  p °.  The  equation 
relating  the  increase  in  volume  to  the  energy  liberated  in  the  reaction 
can  be  readily  derived  by  combining  the  constant  enthalpy  condition  with 
Eqs.  (5)  and  (6).  The  initial  condition  for  propagation  of  the  flame 
along  the  tube  is  (DF/Dt)  > 0.  Other  initial  conditions  follow  from 
the  equations  m = p^(F-u^)  = p^CF-u^.),  and  p - p^.  = m(u.-  u^.) , the 
identities  for  the  pressure  and  particle  velocity  derivative,  and  the 
flow  equations,  Fqs.  (26)  and  (27).  Differentiating  m and  imposing  the 
initial  and  rear-boundary  conditons  leads  to  the  equations 


o 

Dm 

Dt 


o ,DF  °Ui  ° 

Pi  Dt  Dt 


o DF° 
Pf  Dt 


(29) 
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and  rearranging  Eq.  (29)  gives  the  equation 


, .o  DF°_  o 

(pi  Pf}  Dt  “ Pi  Dt 


(30) 


Because  p^  > p^.  , and  (Du./Dt)°  = (du^/dt)°,  it  follows  from  Eq . (30) 
that  (Du^/Dt)  = (du^/dt)  > 0 and  that  the  particle  velocity  In  front 
of  the  flame  starts  to  increase  as  the  flame  starts  to  propagate  along 
the  tube.  That  is,  the  deflagration  process  produces  a compression  wave 
ahead  of  the  flame  front.  Differentiating  the  momentum  equation  and 
imposing  the  initial  and  rear-boundary  conditons  leads  to  the  equations 


o ^o.o,o 
Dpj  Dpf  dpt  dpf 

Dt*  = Dt~  = dt“  = dt~  ( * 

which  shows  that  the  initial  time  rate  of  change  of  pressure  is  the 
same  in  burnt  and  unburnt  materials.  The  assumption  that  material  in 
front  of  the  flame  is  compressed  along  an  isentrope  is  sufficient  to 
determine  other  initial  conditions.  Comhination  of  the  isentropic 
condition  dp^  = (pc)^  du^  with  Eqs.  (30)  and  (31)  gives  the  equation 


, o o 

dpt  dpf 

dt  dt 


o. 

(P, 


.o  DF 
Pf}  Dt 


o 


(32) 


which  shows  that  the  pressure  starts  to  increase  behind  and  in  front  of 
tlie  flame  as  the  flame  starts  to  propagate  along  the  tube.  The  correspond- 
ing equations  for  the  particle  velocity  gradients  obtained  by  combining 

2 

Eq.  (32)  with  the  energy  equation  dp/dt  = - pc  <)u/<)x  show  that  tlie 
initial  values  of  the  particle  velocity  gradients  behind  and  in  front 
of  the  deflagration  are  negative. 

Initial  attempts  to  use  group  theoretical  methods  to  construct 
solutions  for  the  development  of  the  flow  produced  by  the  propagation 
of  a deflagration  wave  in  a closed  tube  have  been  unsuccessful. 
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KXPER I MENTAL  STUDIES 

The  experimental  studies  have  been  delayed  by  problems  of  finding 
a source  of  propellant  for  the  Lagrange  gage  experiments.  Protracted 
negotiations  to  obtain  HMX-based  propellant  from  Lawrence  Livermore 
Laboratory  were  finally  unsuccessful  for  the  following  reasons: 

• Propellant  was  in  short  supply 

» The  machining  costs  and  times  to  meet  tolerances 
in  the  dimensions  of  propellant  samples  needed 
to  construct  the  DDT  gage  assemblies  were 
prohibitive. 

We  subsequently  found  a source  of  HMX  propellant  aL  Edwards  Air  Force 
Base  (EAFB),  and  personnel  there  agreed  to  provide  the  propellant  required 
for  tills  program.  After  discussions  with  EAFB  personnel,  we.  concluded  that 
the  usual  method  of  constructing  targets  for  Lagrange,  gage  experiments 
from  machined  pieces  was  not  practicable  for  HMX-based  propellants.  Hence, 
a target  assembly  based  on  casting  rather  than  machining  the  propellant 
was  designed  to  eliminate  the  complex  machining  and  grooving  operations 
conventionally  used  in  constructing  targets  for  Lagrange  gage  experiments. 
Target  assemblies  consisting  of  acrylic  tubes  containing  the  stress  gages 
and  ionization  pins  are  being  constructed  and  will  be  sent  to  EAFB  to  be 
filled  with  propellant  as  soon  as  they  are  completed.  Thick-wall  steel 
confinement  tubes  to  contain  these  target  assemblies  are  being  constructed 
at  SRI.  We  expect  to  assemble  and  perform  the  initial  1)1)1  experiments 
shortly  after  we  receive  the  cast  charges  from  EAFB. 

FUTURE  WORK 

Theoretical  and  experimental  studies  will  be  continued  during  the 
second  year  of  the  contract  to  model  and  establish  conditions  for  the  onset 
and  occurrence  of  DDT  in  HMX-based  propellants.  Solutions  for  the  unsteady 
flow  produced  by  a deflagration  will  be  sought  to  provide  a basis  for 
interpreting  the  results  of  the  DDT  experiments.  The  pressure  histories 
recorded  in  the  DDT  experiments  will  be  incorporated  into  the  treatment 
of  the  unsteady  flame  to  generate  a satisfactory  model  for  DDT. 


